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Abstract 


This  memorandum  outlines  a  method  for  estimating  the  near-bottom  loads  on  a 
submarine  submerged  in  deep  water  and  for  the  related  problem  of  predicting 
submarine  squat  in  deep  or  shallow  water.  The  method  is  fully  three-dimensional,  but 
has  some  limitations  associated  with  potential  flow  and  linearization  of  the  free 
surface.  The  chief  difficulty  for  this  study  was  the  dearth  of  suitable  data  available  for 
validation.  Nevertheless,  the  results  are  encouraging.  Overall,  it  appears  that  the 
sinkage  component  of  squat  is  predicted  with  greater  accuracy  than  is  trim.  Further 
work  should  include  the  acquisition  of  more,  reliable,  model  data  and  the  investigation 
of  scale  effects. 


Resume 


Le  present  document  donne  un  bref  apergu  d’une  methode  d’estimation  des  charges 
pres  du  fond  d’un  sous- mar  in  en  plongee  en  eaux  profondes  et  des  problemes  connexes 
de  prediction  de  Paccroupissement  d’un  sous-marin  en  eaux  profondes  ou  peu 
profondes.  La  methode  est  entierement  tridimensionnelle,  mais  comporte  certaines 
limites  relatives  a  l’ecoulement  potentiel  et  a  la  linearisation  de  la  surface  libre. 

Malgre  la  rarete  des  donnees  appropriees  disponibles  pour  validation,  principale 
difficulty  rencontree  lors  de  l’etude,  les  resultats  sont  encourageants.  Dans  l’ensemble, 
le  composant  enfoncement  de  Paccroupissement  semble  etre  plus  facile  a  predire  que 
Passiette.  Les  travaux  a  venir  devraient  comprendre  P acquisition  de  plus  de  donnees 
de  modelisation  fiables  et  P analyse  des  effets  de  changement  d’echelle. 


DRDC  Atlantic  TM  2003-078 


i 


This  page  intentionally  left  blank. 


DRDC  Atlantic  TM  2003-078 


Executive  Summary 


Introduction 

Prompted  by  interest  in  the  hazards  associated  with  submarine  operation  in  confined 
water,  an  exploratory  study  was  made  of  a  potential  flow-based  method  of  predicting 
near-bottom  loads  and  squat.  This  memorandum  outlines  the  method  and  some 
validation  with  experimental  data. 


Significance 

A  submarine  operating  near  the  sea  bottom  experiences  a  downforce  and  pitching 
moment  that  complicates  maintaining  bottom  clearance.  Squat,  the  increase  in  draft 
when  underway  on  the  surface,  is  increased  by  these  loads  when  in  shallow  water. 
Understanding  and  predicting  these  effects  is  important  to  the  safe  operation  of  the 
Victoria  class. 


Principal  Results 

The  estimation  method  described  here  predicts  near-bottom  loads  and  squat 
reasonably  well  compared  with  the  few  suitable  validation  data  that  were  found. 
However,  the  trim  component  of  squat  appears  to  be  difficult  for  this  and  other 
methods  to  predict  with  accuracy,  and  furthermore  is  not  easy  to  measure.  More, 
reliable,  experimental  data  are  needed.  The  effects  of  scaling  model  experiments  to  full 
scale  should  be  investigated  further. 


Further  Investigations 

It  is  planned  to  do  near-bottom  experiments  with  a  model  of  the  Victoria  class  in  2003. 


M.  Mackay,  2003,  Estimation  of  Submarine  Near-Bottom  Hydrodynamic  Loads  and 
Squat,  DRDC  Atlantic  TM  2003-078. 

Defence  R&D  Canada  -  Atlantic. 


mi 


DRDC  Atlantic  TM  2003-078 


Sommaire 


Introduction 

Etant  donne  Pinteret  porte  aux  risques  rattaches  aux  deplacements  des  sous-marins  en 
eaux  restreintes,  une  etude  preliminaire  portant  sur  une  methode  de  prediction,  basee 
sur  l’ecoulement,  des  charges  pres  du  fond  et  de  Paccroupissement  a  ete  rnenee.  Le 
present  document  donne  un  aperqu  de  la  methode  et  de  la  validation  obtenue  a  l’aide 
de  donnees  experiment  ales. 


Importance 

Un  sous-marin  qui  se  deplace  pres  du  fond  de  la  mer  est  soumis  a  une  portance 
negative  et  a  un  moment  de  tangage  qui  nuisent  au  maintien  du  degagement  sous  la 
coque.  En  eaux  peu  profondes,  ces  forces  augmentent  Paccroupissement,  qui  equivaut, 
en  surface,  a  l’augmentation  du  tirant  d’eau  lors  du  deplacement.  II  est  important  de 
comprendre  et  de  predire  l’effet  de  ces  forces  pour  assurer  la  securite  de 
fonctionnement  des  sous-marins  de  classe  Victoria. 


Principaux  resultats 

La  methode  d’estimation  decrite  dans  le  present  document  permet  de  predire  assez 
bien  les  charges  pres  du  fond  et  Paccroupissement,  si  l’on  se  fie  aux  quelques  donnees 
de  validation  pertinentes  disponibles.  Cependant,  Pelement  d’assiette  de 
Paccroupissement  semble  difficile  a  predire  avec  precision,  autant  selon  cette  methode 
que  selon  d’autres,  et  est  difficile  a  mesurer.  II  faudrait  plus  de  donnees  experimentales 
fiables.  Une  analyse  plus  poussee  des  effets  de  changement  d’echelle,  essais  sur 
maquettes  a  essais  en  vraie  grandeur,  doit  etre  faite. 


Travaux  futurs 

Des  experiences  sur  les  charges  pres  du  fond  devraient  etre  effectuees  en  2003  avec  un 
modele  de  sous-marin  de  classe  Victoria. 


M.  Mackay,  2003,  Estimation  of  Submarine  Near-Bottom  Hydrodynamic  Loads  and 
Squat,  DRDC  Atlantic  TM  2003-078. 
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1  Introduction 


This  memorandum  discusses  a  method  for  estimating  near-bottom  hydrodynamic 
loads  on  a  submarine  submerged  in  deep  water  and  for  the  related  problem  of 
predicting  squat  for  a  submarine  on  the  surface  in  deep  or  shallow  water.  The 
investigation  was  first  prompted  by  concerns  for  confined  water  operation  of  the  CF 
Victoria  class.  While  those  concerns  have  largely  been  allayed,  the  work  described 
here  may  suggest  worthwhile  directions  for  further  research  and  development. 

The  near-bottom  downforce  experienced  by  a  submerged  submarine  may  present  a 
manoeuvering  problem  since,  without  a  hydrostatic  restoring  force  and  moment, 
control  surface  deflections  are  required  to  maintain  clearance  from  the  sea  floor.  On 
the  surface,  squat  is  the  combination  of  sinkage  and  trim  of  a  ship  submarine  while 
underway.  It  brings  a  likelihood  of  grounding  in  shallow  water,  where  the  additional 
downforce  due  to  proximity  of  the  bottom  exacerbates  the  problem.  In  both  cases  the 
forces  involved  can  be  calculated  by  integrating  dynamic  pressure  over  the  wetted 
surface  of  the  vessel,  and  squat  can  then  be  obtained  from  basic  hydrostatics.  In 
principal  a  similar  procedure  can  be  applied  to  canal  bank  and  ship-ship  interactions, 
although  accurately  modeling  the  free  surface  becomes  an  additional  concern.  Of  the 
two  topics  considered  here,  squat  has  received  the  most  attention  in  the  literature. 

Using  a  number  of  standard  methods,  George  [1]  has  predicted  squat  for  current  CF 
surface  warships;  he  also  summarises  and  compares  the  methods  employed.  Most  of 
them  are  semiempirical  or  entirely  empirical  and,  nominally  at  least,  restricted  to  a 
limited  range  of  hullform  parameters.  Sierra  et  a  1  [2]  discuss  the  effect  of  variable 
depth,  and  provide  a  comprehensive  bibliography  of  squat  and  other  aspects  of 
operation  in  confined  water.  Tuck  [3]  authored  a  general  analytical  solution  based  on 
potential  flow  about  a  slender  hull  in  shallow  water  (with,  however,  some  limitations 
regarding  slenderness  and  shallowness)  that  provides  the  general  form  of  most  of  the 
semiempirical  methods: 

F2 

°  =  (1) 

where  s  is  squat  at  some  location,  e.g.,  AP,  midships,  or  FP;  FI  is  a  function  of  vessel 
geometry;  and  G  is  a  function  of  the  depth  Froude  number  =  U / \fgh,  where  U  is 
ship  speed  and  h  is  depth  of  the  undisturbed  water.  G  is  finite  as  Fh  0,  and 
typically  equal  to  about  1.0.  Equation  (1)  applies  to  the  region  of  practical  interest  for 
confined  water  operation,  the  subcritical  region  where  F^  <  1.  Tuck’s  method  has 
been  extended  in  a  number  of  papers,  e.g.,  reference  [4],  Gourlay,  a  student  of  Tuck, 
further  built  on  this  work  in  an  extensive  review  of  analytical  methods,  including 
calculation  in  the  transcritical  regime  and  for  non-uniform  depth  [5]. 

Standard  methods  for  predicting  squat  generally  do  not  explicitly  break  the  problem 
down  as  suggested  above,  i.e.,  prediction  of  the  hydrodynamic  forces  followed  by  a 
hydrostatic  calculation,  and  are  not  applicable  to  the  estimation  of  near-bottom 
downforce  on  a  submerged  submarine.  Tuck  and  derivative  methods  are  an  exception, 
but  they  are  restricted  to  shallow  water:  h/T  <  2,  where  T  is  draft. 
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The  procedure  described  here  breaks  the  problem  down  as  required,  with  forces 
obtained  by  solving  for  the  potential  flow  numerically,  rather  than  analytically,  using 
the  classic  boundary  integral  method  of  a  Rankine  source  distribution  on  the 
discretized  hull  [6].  An  outline  of  the  calculations  is  given  in  section  2;  however,  a 
number  of  their  attributes  are  summarized  below: 

•  Vessel  representation  is  fully  three-dimensional  —  slenderness  is  not  an 
assumption.  Speed  is  constant. 

•  Potential  flow  inherently  neglects  boundary  layer  growth  and  separation,  but 
corrections  can  be  made  for  thin  boundary  layers. 

•  Canal  banks  and  other  vessels  are  not  represented. 

•  The  free  surface  (in  calculating  forces  on  the  surfaced  submarine)  is  represented 
as  a  symmetry  plane,  or  rigid  streamline;  this  is  a  low  Froude  number  assumption 
analogous  to  the  double- hull  ship  models  sometimes  used  in  wind  tunnels. 

•  The  sea  bottom  is  represented  by  a  symmetry  plane  in  which  the  vessel  is 
reflected;  any  uniform  depth  of  water,  up  to  infinite  depth,  can  be 
accommodated.  However,  the  scope  of  this  study  did  not  permit  modeling  the 
multiple  reflections  that  occur  with  both  bottom  and  surface  symmetry  present, 
so  they  were  neglected. 

•  The  forces  are  not  recalculated  accounting  for  squat  —  a  couple  of  trials 
suggested  that  this  was  not  necessary.  However,  the  hydrostatics  of  surfaced 
submarines  should  properly  account  for  squat  as  noted  below.  Ideally,  the  whole 
sequence  of  calculations  would  be  done  iteratively. 

While  this  approach  provides  almost  all  the  capability  suggested  by  George  for  a  3-D 
method  [1],  it  is  an  exploratory  study  and  the  various  components  of  the  calculation 
have  not  been  optimized  or  organized  into  an  integrated  package.  Ways  of  mitigating 
some  of  the  limitations  and  assumptions  are  discussed  in  the  following  sections. 

There  are  very  few  relevant  data  in  the  open  literature  for  validating  either 
near-bottom  or  squat  calculations  for  submarines,  and  the  amount  of  validation  done 
for  the  present  study  was  limited.  Some  comparisons  of  the  suggested  method  with 
data  and  with  other  predictions  are  presented  in  section  3  for  the  near-bottom  force, 
and  in  section  4  for  squat. 


2  Calculation  Methods 


2.1  Potential  Flow 

Potential  flow  is  calculated  with  the  PC967  suite  of  programs,  which  is  a  revision  of  the 
first-order  boundary  integral  code  EN967  acquired  some  years  ago  for  general  3-D  flow 
calculations  [7]  and  integrated  into  a  more  extensive  ship  hull  flow  procedure  [8].  The 
latter  has  been  superceded  by  the  HLLFLO  code  [9]. 
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The  vessel  is  discretized  into  quadrilateral  and  triangular  panels  each  having  an 
associated  source  strength  and  a  control  point  at  which  the  boundary  condition  of 
tangential  flow  is  enforced.  The  flow  potential  <f>,  which  is  a  function  of  the  unknown 
source  strengths,  must  satisfy  the  Laplace  equation  in  the  flow  domain,  V2,f>  =  0.  The 
equivalence  between  potential  flow  and  the  discrete  model  is  sketched  in  figure  1;  a 
detailed  discussion  of  boundary  integral  methods  can  be  found  in  Hunt  [6] . 

PC967  treats  the  vessel  geometry  as  a  series  of  components.  These  can  include  lifting 
appendages,  i.e. ,  sternplanes,  as  well  as  the  hull.  However,  since  the  effect  of  a 
boundary  can  be  neglected  for  a  lifting  surface  more  than  one  chord  length  away,  only 
the  hull  was  modeled.  Figure  2  shows,  for  a  notional  hull  section,  how  this  was  done 
for  the  three  configurations  of  interest  here.  To  simplify  the  figure,  the  waterline  is 
shown  at  the  equator  of  the  section;  if  this  is  not  the  case,  the  geometry  may  have  to 
be  modified  as  discussed  below.  The  components  coloured  black  in  the  figure  are  input 
explicitly,  and  the  green  image  components  are  accounted  for  by  making  the  bottom  a 
plane  of  symmetry. 

A  slightly  modified  approach  can  be  taken  with  channel  flows  by  discretizing  not  only 
the  wetted  surface  of  the  vessel,  but  also  the  channel  bottom  and  sides,  and  the  free 
surface.  The  flowfield  being  a  finite  rather  than  infinite  domain  introduces  some 
numerical  complications  including  Kelvin  sources  on  the  free  surface  and  the 
requirement  to  match  conditions  at  the  waterline.  If  undertaken,  it  would  be  most 
effective  to  implement  this  capability  in  HLLFLO  [9]  rather  than  the  older  codes. 

Results  are  presented  here  in  standard  submarine  coordinate  axes  and  notation;  see 
the  Nomenclature  and  reference  [10].  In  particular,  vertical  force  and  displacement  are 
positive  downwards,  and  pitch  angle  is  positive  bow  up  (i.e.,  sinkage  and  trim  by  the 
stern  are  both  positive).  This  coordinate  system  is  different  from  that  used  in  PC967 
and  from  that  generally  used  in  ship  hydrodynamics,  so  some  care  is  required  when 
comparing  with  results  from  other  sources. 


2.2  Hydrostatics 


The  squat  of  a  conventional  ship  hull  in  response  to  a  vertical  force  Z  and  pitching 
moment  M  can  be  estimated  using  linearized  hydrostatics  [11].  Sinkage  A z,  which  is 
implicitly  determined  at  the  LCF  (xcf),  and  trim  A 6  are  given  by 


Az 


Z 

pgAWP 


(2) 


A  6 


M 

pg  II 


(3) 


where  Awp  is  the  waterplane  area  and  Ip  is  the  waterplane  moment  of  inertia;  M  and 
Ip  are  relative  to  the  LCF.  The  principal  assumption  is  that  the  vessel  is  essentially 
wall-sided  at  the  waterline.  This  is  indeed  more  likely  to  be  the  case  for  ships  than  for 
surfaced  submarines.  In  particular,  determination  of  Ad  with  equation  (3)  implies 
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that  the  vessel  is  wall-sided  at  the  fore  and  aft  ends,  and  consequently  the  linearized 
formulation  is  questionable  for  a  number  of  configurations  including  surfaced 
submarines  with  the  fore  and  aft  ends  submerged  and  surface  ships  with  shallow 
transom  sterns. 

From  nonlinear  hydrostatics,  the  sinkage  and  trim  are  solutions  of  the  simultaneous 
integral  equations 

rFP 

Z  +  pg\7  =  pg  A(x,Az,A0)dx  (4) 

Jap 

rFP 

M  —  pg\7xcp  =  pg  /  xA(x,  Az,  A6)dx  (5) 

Jap 

where  V  is  the  zero-speed  immersed  volume  and  A(x,  Az ,  AO)  is  the  immersed  section 
area  distribution  for  sinkage  Az  and  trim  AO.  Equations  (4)  and  (5)  can  be  solved 
straightforwardly  with  standard  numerical  techniques.  For  this  study,  the  hydrostatic 
calculations  were  programmed  with  Mathcad;  example  worksheets  are  given  in 
annexes  A  and  B. 


3  Comparison  of  Near- Bottom  Loads 


3.1  Sphere  Moving  Parallel  to  a  Boundary 


PC967  calculations  were  initially  done  for  an  isolated  sphere,  i.e.,  in  an  infinite  flow 
field,  using  288  and  648  panels  (corner  points  every  15  and  10  degrees  respectively). 
These  configurations  are  shown  in  figure  3.  Because  of  symmetry,  the  surface  pressure 
coefficient  cp  is  a  function  of  the  x  coordinate  alone.  The  classical  doublet  model  [12], 
which  is  analytically  exact  for  potential  flow  about  an  isolated  sphere,  gives 


cp(x )  =  l~\ 


—  a  <  x  <  a 


(6) 


where  a  is  the  radius  of  the  sphere.  Predictions  of  cp  are  in  very  good  agreement  with 
the  analytical  solution,  figure  4.  Use  of  a  higher  panel  density  is  not  justified  by  this 
example. 


The  downforce  coefficient  cz  for  a  doublet  equivalent  to  a  sphere  of  radius  a  moving 
above  and  parallel  to  an  infinite  plane  at  height  h  is  equal  to 

cz  =  1(1)  a<h  (7) 

based  on  frontal  area  tt a2  [12],  While  equation  (7)  is  a  good  approximation  for  real 
spheres  at  moderate  clearances,  it  ceases  to  be  so  at  small  clearances  since  the 
boundary  condition  that  the  surface  of  the  sphere  is  a  streamline  is  no  longer  enforced 
by  the  doublet  model.  Using  a  reflected  image  of  the  sphere,  PC967  maintains 
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boundary  conditions  both  on  its  surface  and  at  the  location  of  the  plane,  resulting  in  a 
coefficient  that  correctly  approaches  infinity  as  h/a  approaches  1.0,  figure  5. 

The  rapid  increase  in  downforce  at  small  clearances  is  further  illustrated  by  plotting 
the  pressure  coefficient  round  a  transverse  cross-section  close  to  the  center  of  the 
sphere,  figure  6.  It  can  be  seen  that  the  upper  surface  of  the  sphere  has  a  very  small 
incremental  effect  on  the  downforce,  even  at  very  small  clearances. 

The  characteristics  of  a  somewhat  submarine-like  hull  can  be  illustrated  by  the 
addition  of  parallel  midbody.  Figure  7  depicts  288-panel  spheres  stretched  by  midbody 
of  length  2a,  where  in  one  case  the  midbody  is  represented  as  one  section  with  high 
aspect  ratio  panels  and  in  the  other  by  eight  sections  so  that  the  midbody  panels  have 
similar  dimensions  to  adjacent  panels  on  the  hemispherical  end-caps.  The  axial 
pressure  distribution  for  these  in  isolation  shows  a  drastic  decrease  of  the  suction 
(negative  cp )  peak  at  maximum  radius  compared  with  the  original  sphere,  figure  8. 
However,  a  level  of  reduced  suction  is  maintained  along  the  midsection,  and  the 
accelerated  flow  that  this  represents  is  a  major  contributor  to  generating  a 
near-bottom  force.  Figure  9  illustrates  how  the  downforce  in  the  proximity  of  a  plane 
is  increased  for  a  stretched  vs.  unstretched  sphere. 


3.2  Submarine  Model  Experiments 

In  reference  [13]  Bystrom  and  Andersson  describe  captive  near-bottom  experiments 
with  a  model  of  the  A17  submarine  (Vastergotland  class).  Data  in  that  paper  were 
sanitized  by  omission  of  the  ordinate  axis  scale  from  the  graphs,  so  this  memorandum 
makes  primarily  qualitative  comparisons.  The  PC967  calculations  were  done  using  a 
representative  hull  (the  overall  dimensions  are  in  proportion,  but  the  tail  is  probably 
fuller  than  the  prototype)  with  30  panels  longitudinally  and  24  around,  figure  10. 

Reference  [13]  plots  the  experimental  loads  against  2 r/c,  where  r  is  maximum  hull 
radius  and  c  is  clearance  between  the  keel  and  the  tank  bottom.  Data  for  level  trim, 
corresponding  to  the  upper  and  lower  parts  of  figure  6  in  reference  [13],  are  shown  in 
figures  11  and  12,  where  the  data  and  PC967  calculations  have  been  adjusted  to  give 
reasonable  overall  agreement.  These  figures  also  include  linear  fit  lines  suggested  in 
reference  [13]  for  modeling  the  near-bottom  effect  in  the  equations  of  motion.  Despite 
scales  not  having  been  optimally  matched,  the  form  of  the  PC967  curves  characterizes 
the  data  better  than  the  linear  relationship,  although  neither  captures  very  well  the 
trend  with  2 r/c  seen  in  figure  12. 

The  same  data  and  PC967  predictions  are  plotted  in  figures  13  and  14,  using  the  same 
relative  ordinate  scales  as  before,  against  h/r ,  where  r  is  maximum  hull  radius.  Since 
h/r  is  equivalent  to  h/a  in  section  3.1,  figure  13  is  analogous  to  figures  5  and  9. 
Presented  in  this  fashion,  comparisons  between  the  data  and  potential  flow  appear 
more  reasonable,  while  the  anomalous  behavior  in  measured  pitching  moment  as 
h  — >  r  is  more  evident.  Following  the  form  of  equation  (7),  the  near-bottom  downforce 
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coefficient  can  be  represented  by 


Z'  oc 


(8) 


where  Z'  is  based  on  l?  [10]  and  L  is  hull  length,  with  a  similar  expression  for  the 
pitching  moment  coefficient  M' .  Since  n  is  unaffected  by  scaling  the  loads,  estimates 
of  it  can  be  derived  from  figures  13  and  14: 


Reference  [13]  data 
PC967  prediction 


Z'  M' 

Figure  13  Figure  14 


2.8 

1.2  to  2.2 

2.3 

3.0 

The  indicated  range  of  n  for  M'  data  was  obtained  by  including  or  omitting  the 
anomalous  pitching  moment  measurements.  Ignoring  the  anomaly,  the  near-bottom 
effect  on  a  submarine  can  be  modeled  by  equation  (8)  and  its  M'  equivalent  with,  as  a 
rough  rule  of  thumb,  2  <  n  <  3.  The  potential  flow  method  outlined  here  should  be 
able  to  estimate  both  n  and  the  constant  of  proportionality  reasonably  well,  but 
experimental  corroboration  is  required  to  accurately  model  a  specific  submarine. 


4  Comparison  of  Squat  Results 


4.1  CPF  in  Deep  and  Shallow  Water 

Since  there  are  very  few  suitable  squat  data  for  submarines,  this  section  presents  a 
discussion  of  some  results  obtained  for  the  Canadian  Patrol  Frigate  (CPF). 

The  ship  was  taken  to  be  in  Deep  Departure  condition  with  FP  draft  of  4.99  m  and 
AP  draft  of  4.95  nr,  corresponding  to  the  model  experiments  described  in 
reference  [14].  The  wetted  surface,  Component  1  in  the  notation  of  figure  2,  and  its 
image,  Component  2,  are  represented  in  PC967  by  20  panels  longitudinally  and 
12  panels  around  the  girth;  i.e.,  their  corner  points  are  defined  at  21  equally  spaced 
stations  in  15  degree  increments.  The  paneling  is  shown  in  figure  15.  Hydrostatic 
calculations  are  illustrated  by  annex  A. 

The  simple  squat  model  represented  by  equation  (1)  clearly  does  not  apply  to  deep 
water  since  it  is  zero  in  the  limit  h  — >  oo.  Deep  water  squat  is  therefore  plotted  against 
speed  or  the  conventional  hull  length  Froude  number  Fr  =  U/y/gL.  Figure  16 
compares  some  predictions  with  scaled-up  model  test  data.  The  notable  feature  of  the 
measurements  is  a  reversal  from  bow  down  trim  at  low  speeds  to  much  larger  bow  up 
trim  at  Fr  ~  0.4  ( U  ~  25  kt),  which  is  not  predicted  by  either  the  linear  or  nonlinear 
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hydrostatic  calculations.  The  figure  also  shows  that  the  linear  and  nonlinear 
predictions  give  about  the  same  sinkage  but  that  trim  is  markedly  less  bow  down  in 
the  nonlinear  prediction,  presumably  due  to  the  lack  of  wall-sidedness  at  the  fore  and 
aft  ends. 


In  the  absence  of  data  for  shallow  water  squat,  some  comparisons  were  made  with  the 
methods  discussed  in  reference  [1].  Unfortunately,  the  key  figures  in  that  reference, 
e.g.,  figure  B-l-1,  are  not  suitable  for  direct  comparison  since  the  trim,  and  whether 
maximum  squat  is  at  the  FP  or  AP,  are  not  noted.  The  effect  of  water  depth  on  the 
longitudinal  distribution  of  downforce  calculated  by  PC967  for  the  CPF  is  illustrated 
in  figure  17.  Squat  predictions  were  compared  with  the  methods  of  Vermeer  [15]  and 
Millward  [16]. 


Vermeer’s  method  is  based  on  that  of  Tuck  [3]: 

a  V  Fh 

AzMid  =  cz-p>  / 

L  2  -  Fl 


A  6  =  C0 


V 

13 


Fl 


(9) 


(10) 


where  Azmid  is  sinkage  at  midships.  He  simplified  Tuck’s  analytical  derivation  of  the 
coefficients  cz  and  cq  by  approximating  the  waterplane  curve  B(x )  and  area  curve 
S(x)  by 

(11) 


B(x)  =  (l  —  x2^j  (l  +  as  x2  +  (3b  x3^j 
S(x)  =  (l  -  x 2)  (l  +  as  x2  +  (3s  z3) 


(12) 


in  which  as,  (3b ,  crs,  and  (3s  are  functions  of  LCF,  longitudinal  center  of  buoyancy 
LCB,  and  waterplane  and  prismatic  coefficients,  Cw  and  Cp  respectively.  For  shallow 
water  of  infinite  extent,  Vermeer  obtained: 


cz  =  „ — -p—  [32  —  40(Cp  +  Cw)  +  75  CpCw  —  980  CpCwCbCf]  (13) 

07TC  W'-'P 


ce 


-7 

I87 rCwCpkyy 


(20  Cp  —  45  CpCw)Cp  +  (24  Cw  ~  39  C2y)Cp 


(14) 


where  kw  is  longitudinal  radius  of  gyration  of  the  waterplane  and  Cf  and  Cp  are  the 
LCF  and  LCB  nondimensionalized  by  L.  He  gives  similar  expressions  for  sinkage  and 
trim  in  a  finite  channel  [15].  Vermeer’s  approximations  to  the  CPF  waterplane  and 
area  curves  are  shown  in  figure  18;  it  can  be  seen  that  equations  (11)  and  (12)  assume 
a  cruiser,  rather  than  transom,  stern. 


From  a  large  number  of  towing  tank  tests  with  commercial  ship  models  in  shallow 
water  of  effectively  infinite  extent,  Millward  obtained  empirical  expressions  [16]  for 
squat  at  midships  (smid)  and  at  the  bow  (show)  hr  the  form  of  equation  1: 

SMID  -  (12.22  CBf-  0.46)  (15) 
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SBOw  =  (15.0  CB  §  -  0.55)  r35L_  (16) 

where  B  is  the  beam  at  midships.  The  range  of  validity  of  these  expressions  is  stated 
to  be  1.25  <  h/T  <  6. 

Figures  19  and  20  compare  linear  estimates  of  AP  and  FP  squat  from  the  present 
method  with  calculations  from  Vermeer’s  and  Millward’s  methods.  There  is  clearly  a 
considerable  difference  between  all  three.  Vermeer  gives  somewhat  less  sinkage,  and 
significant  bow  up  trim,  compared  with  the  modest  bow  down  trim  predicted  by 
Millward;  this  is  consistent  with  George’s  figure  3-6  [1].  Since  the  present  method 
accounts  for  depth  and  speed  separately,  not  combined  into  the  single  parameter  F^.  it 
generates  a  family  of  curves  when  squat  is  plotted  against  the  depth  Froude  number, 
as  seen  in  these  figures.  The  “fit  to  linear  estimate”  curves  are  locally  weighted 
regressions  of  both  the  h  =  8  m  and  12  nr  estimates.  These  curves  indicate  a  predicted 
bow  down  trim  similar  to  Millward  with  somewhat  less  sinkage  —  more  like  that 
suggested  by  Vermeer.  Without  data  for  shallow  water  squat  it  is  not  possible  to 
resolve  the  discrepancies  between  these  predictions. 


4.2  Submarine  in  Shallow  Water 

In  reference  [17],  Eddison  and  Fryer  describe  a  study  of  surface  manoeuvering  in 
confined  water  for  the  SSBN05  (Vanguard)  class  submarine.  The  study  included 
model  scale  measurements  of  squat  in  shallow  water. 

The  draft  of  a  surfaced  submarine  is  typically  at  least  80  percent  of  pressure  hull 
diameter,  as  shown  in  the  left-hand  sketch  of  figure  21,  rather  than  50  percent  as 
suggested  by  figure  2.  Paneling  the  hull  and  its  reflected  image  to  the  waterline  causes 
difficulty  for  the  potential  flow  calculation;  errors  arising  from  proximity  of  hull  and 
image  panels  at  the  waterline,  and  from  an  unrealistic  local  surface  boundary 
condition,  may  become  large  relative  to  the  contribution  of  the  upper  hull  as  the  draft 
increases.  To  mitigate  these  problems,  in  potential  flow  calculations  the  lower  hull  was 
joined  to  its  image  by  vertical  panels  at  the  maximum  beam,  as  shown  in  the  center 
and  right-hand  sketches  of  figure  21.  Although  this  effectively  determines  downforce 
and  pitching  moment  from  just  the  lower  hull  and  image  contributions  to  it,  the 
approximation  is  considered  acceptable. 

Figure  22  shows  paneling  arrangements  for  an  axisymmetric  hull  representing  SSBN05 
general  proportions,  and  the  corresponding  modified  hull  and  image  for  surfaced 
downforce  and  pitching  moment  estimation.  Notwithstanding  the  geometric 
modification  required  for  PC967,  hydrostatic  calculations  were  done  with  the  actual 
geometry,  see  annex  B. 

The  bare  hull  geometry  is  normally  sufficient  for  near-bottom  or  squat  force 
calculations.  The  effect  of  a  including  a  keel  was  investigated  for  this  hull  by 
calculating  the  deep  water  near-bottom  force  for  the  bare  hull  and  with  aim  deep, 

2  m  wide  keel  added  along  the  parallel  midbody,  about  57.5  percent  of  total  hull 
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length.  Height  above  the  bottom,  h ,  was  set  at  8.4  m,  giving  a  zero-trim  clearance  of 
2  m  for  the  bare  hull  and  1  m  for  the  keel.  The  resulting  downforce  distributions  are 
shown  on  figure  23.  The  difference  in  net  downforce  is  less  than  four  percent,  while  the 
keel  brings  the  peaks  of  the  distribution  a  little  closer  to  midships,  i.e.,  towards  where 
it  begins  and  ends.  These  differences  do  not  justify  modeling  a  keel  even  at  very  small 
bottom  clearances.  It  was  noted  previously  that  it  is  not  necessary  to  model  cruciform 
configuration  sternplanes,  which  are  about  one  chordlength  above  the  bottom  at 
minimum  clearance;  studies  of  wings  in  ground  effect,  e.g.,  reference  [18],  show  that 
the  interaction  at  this  distance  is  negligible.  Interaction  with  lower  X-rudders,  on  the 
other  hand,  might  not  be  negligible. 

The  shallow  water  squat  data  in  reference  [17]  are  reproduced  on  figure  24.  There  is  a 
considerable  discrepancy  between  data  for  1/40  and  1/26  scale  models.  Eddison  and 
Fryer  concluded  that  the  latter  incorporated  excessive  blockage  effects,  so  they  are  not 
included  in  the  comparisons  below.  It  was  observed  that  the  models  trimmed  by  the 
bow;  however,  the  “estimated  sinkage  at  CG”  curve  on  the  figure  indicates  that  the 
trim  angle  was  quite  small.  Reference  [17]  gives  little  detail  on  experimental 
conditions.  In  the  comparisons  presented  here,  squat  estimates  were  made  for  the 
submarine  initially  in  nominal  deep  draft  condition  with  a  trim  of  0.7  degrees  by  the 
stern. 

Although  a  surfaced  submarine  differs  quite  considerably  from  the  hullforms  on  which 
the  Vermeer  and  Millward  methods  are  based,  it  is  of  interest  to  first  compare 
predictions  from  those  methods  with  the  model  data.  Vermeer’s  approximations  to  the 
waterplane  and  area  curves,  figure  25,  are  reasonable  except  near  the  stern,  where,  as 
the  waterplane  curve  shows,  the  hull  sections  were  fully  submerged.  Comparison  with 
data  for  stern  sinkage  is  shown  in  figure  26.  Forward  sinkage  was  predicted  to  be 
significantly  higher  than  measured,  and,  as  the  mean  sinkage  approximates  the 
estimate  for  the  CG  on  figure  24,  it  appears  that  overprediction  of  trim  by  the  bow  is 
largely  responsible  for  the  discrepancies  in  this  figure. 

Potential  flow  calculations  were  done  at  model  scale  and  full  scale,  modifying  the  hull 
geometry  in  each  case  by  a  boundary  layer  displacement  thickness  distribution  derived 
for  axial  flow  on  the  axisymmetric  hull  [19].  Displacement  thickness  8*  is  compared  for 
model  (1/40)  scale  at  various  full  scale  speeds  and  full  scale  at  2  kt  in  figure  27. 

It  was  noted  in  the  previous  section  that  squat  should  be  represented  by  a  family  of 
curves  when  plotted  against  depth  Froude  number.  The  predictions  in  figure  28 
represent  calculations  done  for  h  =  10,  11,  12,  13,  14,  and  15  m  at  2,  4,  6,  8,  10,  and 
12  kt  full  scale,  and  scaled  up  to  6,  8,  10,  and  12  kt  for  a  1/40  model.  Model  scale 
predictions  are  in  reasonable  agreement  with  the  data,  with  the  full  scale  stern  sinkage 
predicted  to  be  20  to  25  percent  lower. 

Notwithstanding  the  similarity  of  model  scale  and  full  scale  stern  sinkage  predictions 
in  figure  28,  forward  sinkage  agrees  less  well.  This  is  because  trim,  although  quite 
small,  is  predicted  to  be  different  in  the  two  cases,  figures  29  and  30.  In  figure  30,  the 
model  is  predicted  to  trim  by  the  stern  while  full  scale  trims  by  the  bow  at  low  speed 
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and  at  shallow  depths,  otherwise  tending  to  trim  by  the  stern.  Further  corroboration 
is  required  to  establish  whether  this  is  a  genuine  scaling  effect  or  an  artifact  of  the 
calculation  method. 


5  Concluding  Remarks 


This  memorandum  records  a  brief  study  of  estimating  submarine  near-bottom  loads 
and  surfaced  squat  by  combining  potential  flow  hydrodynamic  forces  with  hydrostatic 
calculations.  Squat  estimates  were  compared  with  the  standard  methods  (for  surface 
ships)  of  Vermeer  and  Millward.  The  scarcity  of  suitable  experimental  data  for 
submarines  presents  a  problem  for  adequately  validating  both  near-bottom  and  squat 
calculations. 

Boundary  integral  calculation  of  potential  flow  near-bottom  force  was  validated  for  a 
sphere  moving  parallel  to  the  bottom  and,  with  respect  to  the  trend  with  clearance 
under  the  hull,  for  experiments  with  a  submarine  model,  although  absolute  values  for 
comparison  were  not  available  in  the  latter  case.  DRDC  is  planning  near-bottom 
model  experiments  that  will  enable  more  detailed  comparisons  to  be  made.  The 
principal  test  parameters  will  be  clearance  and  trim  angle;  secondary  parameters  may 
include  sternplane  deflections  and  drift  angle. 

Squat  estimates  were  compared  with  model  data  for  the  CPF  in  deep  water  and  for  a 
submarine  in  shallow  water.  Results  obtained  with  the  present  method  and  other 
methods  discussed  here  suggest  that  sinkage  can  in  general  be  predicted  at  least  more 
consistently,  if  not  more  accurately,  than  can  trim.  Nevertheless,  accurate 
measurement  of  these  quantities,  particularly  trim,  is  not  easy  even  in  well-conducted 
model  experiments,  which  adds  another  element  of  uncertainty  to  comparisons  and 
validation.  Hydrostatic  calculations  for  the  submarine  demonstrated  scaling  effects  in 
squat;  this  should  be  investigated  in  future  investigations  since  standard  estimation 
methods  rely  heavily  on  model  scale  experimental  data. 
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BOUNDARY  INTEGRAL  METHOD 


V  =  V4> 

v2<j>=  o 

POTENTIAL  FLOW 


Y 

= 

f  °'i  °j 

v, 

<T\ 

v2 

= 

°ij 

% 

DISCRETE  MODEL 

Figure  1.  Equivalence  of  potential  flow  and  the  discrete  model,  from  [8].  Freestream 
flow,  equivalent  to  U  in  the  text,  is  denoted  V^;  V  and  Vi  are  local  velocity;  are 
geometrically-determined  influence  coefficients;  and  Uj  are  the  source  strengths. 


Submerged,  Deep  Water  Surfaced,  Deep  Water  Surfaced,  Shallow  Water 


Wetted  Surface 
Primary  Image 
Secondary  Image 


Comp.  1' 


, - Comp.  1' 


Comp.  2' 


Figure  2.  PC967  geometry  setup  for  three  configurations:  submerged,  deep  water; 

surfaced,  deep  water;  and  surfaced,  shallow  water.  For  each,  the  calculation  sums 
the  contributions  of  all  components  to  the  force  on  Component  1. 
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Figure  3.  Sphere  paneling:  left,  288  panels;  right,  648  panels.  The  x  axis  is  horizontal. 


Figure  4.  Axial  pressure  distribution  on  an  isolated  sphere. 
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Figure  5.  Force  coefficient  on  a  sphere  moving  parallel  to  a  plane. 


Figure  6.  Pressure  coefficient  round  the  near-midsection  of  a  648-panel  sphere 

moving  parallel  to  a  plane. 
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Figure  7.  Stretched  spheres  (288+  panels):  left,  “1  panel”  stretch; 
right,  “8  panel”  stretch.  The  x  axis  is  horizontal. 
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Figure  8.  Axial  pressure  distribution  on  isolated  stretched  spheres. 
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0.6 

0.5 

0.4 

J'1  0.3 

0.2 

0.1 

0.0 

Figure  9.  Force  coefficient  on  a  sphere  and  stretched  sphere  moving  parallel  to  a  plane. 


288  panel  sphere 
8  mid  section  stretch 


Figure  10.  Representative  hull  for  comparison  with  near-bottom  model  experiments  [13]. 
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Figure  11.  Submarine  near  bottom:  downforce  plotted  against  2 r/c, 
as  in  reference  [13],  figure  6. 


Figure  12.  Submarine  near  bottom:  pitching  moment  plotted  against  2 r/c, 

as  in  reference  [13],  figure  6. 
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h/r 


Figure  13. 


Submarine  near  bottom:  downforce  plotted  against  h/r. 


Figure  14.  Submarine  near  bottom:  pitching  moment  plotted  against  h/r. 
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Figure  15.  Paneling  for  the  CPF  hull  and  its  image:  top,  side  elevation; 
below,  plan  view;  left,  section  plan;  right,  isometric  view. 
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Figure  16.  CPF  deep  water  squat.  The  lines  are  predictions,  the  symbols 
are  data  from  reference  [14]. 
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-10 


Figure  17.  CPF:  longitudinal  distribution  of  downforce  calculated  by  PC967  for  h  =  oo, 
12  m,  and  8  m.  Here,  and  following,  station  0  =  FP  and  station  20  =  AP. 


Figure  18.  Nondimensional  CPF  beam  and  area  compared  with  Vermeer's 
approximations  (equations  (11)  and  (12)). 
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Squat,  m,  at  FP  Squat,  m,  at  AP 


0.0 


0.1 


0.2 


0.3 


0.4 


• 

Linear  est.,  h  =  8  m 

o 

Linear  est.,  h  =  12  m 

Fit  to  linear  est. 

Vermeer 

Milward 

\  \ 


_ A _ 


—I _ 


0.0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 
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Figure  19.  CPF:  shallow  water  AP  squat  estimates. 


Figure  20.  CPF:  shallow  water  FP  squat  estimates. 
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Actual 


Deep  Water 


Shallow  Water 


-  Wetted  Surface 

.  Non-wetted  Surface 

-  Primary  Image 

.  Secondary  Image 


Bottom 

Comp.  1' 

\ 


Comp.  2' 


Figure  21.  Actual  geometry  for  a  surfaced  submarine  hull,  and  PC967  approximations 
for  deep  and  shallow  water.  The  calculation  sums  the  contributions  of 
all  components  to  the  force  on  Component  1. 


Figure  22.  Submarine  hull  paneling  for  comparison  with  squat  experiments: 
left,  submerged,  24  x  40  panels;  right,  surfaced  (with  image),  36  x  40  panels. 
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Figure  23.  Longitudinal  distribution  of  downforce  estimated  by  PC967  for 
SSBN05  hull  with  and  without  a  keel  in  deep  water  at  h  =  8.4  m. 


Figure  24.  Submarine  squat:  sinkage  data  for  SSBN05,  from  reference  [17]. 
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Figure  25.  Nondimensional  SSBN05  beam  and  area,  in  nominal  surfaced  deep  draft 
condition,  compared  with  Vermeer's  approximations  (equations  (11)  and  (12)). 


Figure  26.  Submarine  squat:  comparison  of  SSBN05  model  data  [17]  for  stern  sinkage 
with  estimates  from  the  Vermeer  and  Millward  methods. 
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Figure  27. 


Estimated  nondimensional  boundary  layer  displacement  thickness  for 
SSBN05  hull  in  axial  flow  (deep  water). 


Figure  28.  Submarine  squat:  comparison  of  SSBN05  model  data  [17]  for  stern  sinkage 
with  estimates  at  model  and  full  scale  from  the  present  method. 
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Figure  29.  Submarine  squat:  predicted  sinkage  at  LCF  for  SSBN05. 


Figure  30.  Submarine  squat:  predicted  trim  for  SSBN05. 
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Annex  A.  Hydrostatic  Calculation  Worksheet  for  Surface  Ship 


With  minor  modifications,  the  procedure  shown  on  the  next  five  pages  could  be  used 
for  any  surface  ship.  The  worksheet  is  programmed  in  Mathcad  2000 
( www . mathcad . c  om) . 

Hull  geometry  to  the  sheer  line  is  entered  into  the  input  tables  YS  and  ZS  for  21 
stations,  starting  at  the  AP.  ZS  is  measured  up  from  the  keel  line.  Right-click  within 
each  table  and  select  “Import”  to  read  an  ASCII  hie  of  values.  Normal  force  and 
pitching  moment  from  PC967,  which  are  entered  at  the  top  of  the  third  page  of  the 
worksheet,  must  be  in  standard  submarine  coordinates  (Z  downwards,  M  bow  up). 

Linear  squat  estimates  are  made  with  the  calculation  of  two  intermediate  quantities: 
TPI  (Tons  Per  Inch  immersion  -  the  name  is  a  historical  anachronism)  and  MCT 
(Moment  to  Cause  Trim)  [11],  The  sinkage  and  trim  are  Z/ TPI  and  M/ MCT 
respectively. 

In  the  iterated  estimates,  note  that  section  area  is  obtained  by  integrating  the 
interpolated  section  offsets.  Consequently,  the  right-hand  sides  of  equations  (4)  and 
(5),  which  are  on  the  last  page  of  the  worksheet,  are  double  integrals  with  internal 
computational  overhead  for  interpolation  and  several  function  evaluations.  This  makes 
solving  for  Ac  and  Ad  very  slow,  and  relatively  few  of  these  evaluations  were  done  for 
the  present  study.  The  process  would  be  considerably  speeded  up  by  generating  a  set 
of  Bonjean  (section  area)  curves  prior  to  the  calculation. 
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HALHULL  STAT2.MCD 


Section  Geometry: 

YS  := _ 


0 

1 

2 

3 

4 

5 

6 

0 

0 

0 

0 

0 

0 

0 

0 

1 

2000 

1000 

148 

178 

341 

1000 

1000 

2 

3500 

2000 

1000 

1000 

759 

1183 

2000 

3 

5000 

3000 

2000 

2000 

1000 

2000 

3000 

4 

5434 

3871 

3000 

3000 

1716 

3000 

3437 

5 

5549 

4000 

4000 

3293 

2000 

3085 

4000 

6 

5630 

5000 

5000 

4000 

3000 

4000 

4896 

7 

5767 

5807 

5131 

5000 

4000 

4787 

5000 

8 

5924 

5913 

6000 

5568 

5000 

5000 

5777 

9 

6000 

6000 

6117 

5892 

5376 

6000 

6000 

23  := 


0 

1 

2 

3 

4 

5 

6 

0 

4630 

3670 

2970 

1880 

517 

0 

0 

1 

4631 

3739 

3000 

2000 

1000 

897 

345 

2 

4632 

3811 

3153 

2453 

1500 

1000 

647 

3 

4633 

3894 

3318 

2750 

1682 

1281 

895 

4 

5000 

4000 

3488 

2944 

2000 

1483 

1000 

5 

5220 

4020 

3675 

3000 

2075 

1500 

1155 

6 

5500 

4286 

3945 

3152 

2260 

1720 

1500 

7 

6000 

5000 

4000 

3461 

2465 

2000 

1549 

8 

7000 

5220 

4677 

3750 

2803 

2098 

2000 

9 

7486 

5461 

5000 

4000 

3000 

2793 

2176 

Some  definitions: 


0..20 

YS 

23 

Scale  section  data  to  m 

1000 

1000 

124.5 

p  :=  1025 

g:=  9.81 

Input  initial  draft  (to  hull  keel  line)  at  AP  (first)  and  FP  (second): 


Ts  := 


0 

0 

4.950 

1 

4.990 

AP 

TA  :=  Tsj. 

FP 

TF  :=  Tsj 
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Function  Statements 


T(J,Ta,Tf)  :=  Ta  +  —  (Tf  -  Ta) 


Draft 


Y(Z,J) 


:=  linterpt 


zs<j\ys<j\z) 


ywpj  :=  2  Y(T(j,TA,TF) 


ywpj  :=  if|ywpj  >  0,ywpj,0j 


xh.  :=  — (j  -  10) 

J  20  U  } 


y(z) 

Waterplane  breadth 

ywp20  :=  o 

Axial  coordinate  wrt  midships,  x fwd 


20 


ywp 


10  - 


-80 


Ywp(x)  :=  linteip(xh,ywp,x) 


^xh- 


AwP  := 


■20 


Ywp(x)  dx 


xho 

TPI  :=  AwP'P'g 
rxhoo 


Awp  =  1437.7788 


TPI  =  1.4457  x  10 


x  Ywp(x)  dx 


LCF  := 


xho 


Awp 


LCF  =  -8.0372 


Waterplane  area 

TPI  in  newtons  per  metre 


wrt  midships  (negative  -  aft) 


rxh20 

IL- 

J 

(x  -  LCF)2  Ywp(x)  dx 
xho 

lL=  1.3483  x  10 

Il  p  g  n 

MCTa  - - 

180 

MCT(j  =  2.3662  x  108 

MCT  in  Newton-m  per  degree 

iL'Pg 

MCTm  := 

L 

MCTm=  1.0889  x  108 

MCT  in  Newton-m  per  metre 

80 
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Input  normal  (down)  force  and  pitching  moment  about  midships 


0 

0 

3.58-106 

1 

-1.57-107 

Z 

M 


Shift  moment  to  LCF: 


Linear  estimates  of.. 


Z  :=  zm^ 
M  :=  zm, 


Mlcf  :=  M  +  Z-LCF 


6 


Z  =  3.5800  x  10 


M  = -1.5682  x  10 


Mlcf  =  -4.4455  x  10 


Sinkage: 


Trim: 


Az  := 


A©  := 


TPI 

Mlcf 

MCTd 


Az  =  0.24763 


A@  =  -0.18788 


mLCF 

MCTm 


metres,  at  LCF 

degrees 


=  -0.40825  metres,  positive  by  stern 


Volumetric  Properties: 


area.  :=  if 
J 


ZSQ  .  <T(j,TA,TF),2- 


rT(j,TA,TF)  s' 

Y(z,j)  dz,0 


Area(x)  :=  linterp(xh,area,x) 


Vol  := 


A, 


3*0 


Area(x)  dx 


ZS, 


o,j 


(Note,  trapezoidal  integration  will  underestimate 
volumetric  properties  by  a  little  bit.) 


Vol-p 

1000 


:  4630.6155 


Salt  water  displacement,  t. 
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Zmit ;=  P'g  Vol 


Zjjut  =  4.5426  x  10 


Initial  Z  (hydrostatic  force) 


/•xho. 


MomV  := 


x-Area(x)  dx 


xho 

Mmit :=  p  g-MomV 

.  __  MomV 
LCB  := 


Vol 


Mjnit  =  -l  0861  x  10s  Initial  M  (hydrostatic  moment) 


LCB  =  -2.3910 


(relative  to  midships) 


Iterated  Estimates: 


Definitions: 


Zt  :=  Zjait  +  z  Initial  +  Squat  Force 
and  moment 

:=  Miirit  +  M 


TAn(dz,d©)  :=  TA  +  dz  +  ^LCF  +  ^  tanfde— 


TFn(dz,de)  :=  TF  +  dz  +  ^LCF  -  ^  j-tan|^d9"^ 


Note: 


Zinit 


■  0.0788 


M 

Mmit 


=  0.1444 


At  AP:  TAn(Az,Ae)  -  TA  =  0.0699  m 
At  FP:  TFn(Az,Ae)  -  TF  =  0.4781  m 


(Linear  squat) 


Tn(j,dz,de)  :=  TAn(dz,de)  +  —  (TFn( dz.de)  -  TAn(  dz.de)) 


20 


are  an(j,  dz.de)  :=  if 
Arean(x,  dz.de)  := 


■.T^J.dz.de) 


ZSQ  J<Tn(j,dz,de),2-  Y(c,j)d^,0 

,  ’  Jzs0iJ 


for  k  e  0  . 20 
A^  <-  are  an(k,  dz.de) 

linterp(xh.A.x) 


6z  :=  Az  69  :=  AO 
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Solve  Loop 


Given 


Zt=  p  g 


,-xh. 


l20 


Arean(x,6z,6©)  dx 


Mj  =  P  g 


xho 

»xh20 

xho 


xArean(x,&z,S6)  dx 


f  Az 
<A9 


Find(Sz,Se) 


(may  be  disabled) 


Az  =  0.2457  m  (at  LCF) 

A9  =  -0.0530  degrees 


Iterated  Squat 

At  AP:  TAn(Az,A9)  -  TA  =  0.1955 

At  FP:  TFn(Az,Ae)  -  TF  =  0.3108 
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Annex  B.  Hydrostatic  Calculation  Worksheet  for  Submarine 


With  minor  modifications,  the  procedure  shown  on  the  next  four  pages  could  be  used 
for  any  surfaced  submarine  or  vehicle  for  which  the  hull  is  axisymmetrical. 

Longitudinal  coordinate  x,  positive  aft,  and  hull  diameter  are  entered  into  the  input 
tables  hp  for  arbitrary  stations,  starting  at  the  FP.  Right-click  within  each  table  and 
select  “Import”  to  read  an  ASCII  hie  of  values;  a  modified  PC967A  input  hie  was  used 
in  this  instance.  Normal  force  and  pitching  moment  from  PC967,  which  are  entered 
near  the  top  of  the  third  page  of  the  worksheet,  must  be  in  standard  submarine 
coordinates  (Z  downwards,  M  bow  up). 

Linear  squat  estimates  are  made  as  in  annex  A.  For  the  iterated  estimates,  section 
area  is  an  analytical  function  of  local  sinkage.  The  right-hand  sides  of  equations  (4) 
and  (5)  are  therefore  single  integrals,  and  solving  for  /S.z  and  A0  is  much  faster  than 
for  the  worksheet  in  annex  A. 
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AXHULL_STAT_2.MCD 

Input  Table  for  x  and  diameter,  x  going  aft  (can  use  modified  PC967A.INP  file  ): 


0 

1 

0 

-74.9 

0 

1 

-71.35 

7.3 

2 

-67.6 

8.9 

3 

-63.85 

10.2 

4 

-60.1 

11.4 

5 

-56.35 

12.5 

(right-click  a  cell  and 
"import''  to  change  file) 

n  :=  0..  N 


xh  :=  -hp 


N-n,0 


dia  :=  hp,,  , 
n  rN-n,l 


L  :=  xh. 


N  ^ 


N  :=  len  gth(hp<0>)  -  1 
N  =  40.0000 

Assign  columns;  sort 
(N  fwd,  0  aft)  ;setxfwd 


D  :=  max(dia) 
L  =  149.8000 


D  =  12.8000 


Input  Table  for  initial  draft  (to  hull  keel  line)  at  AP  (first)  and  FP  (second): 


Ts  := 


0 

0 

9.915 

1 

8.085 

AP 

TA  :=  Tsj. 

FP 

TF  :=  Tsj 

V=  TA-Q^-xho)- J 

ywp^  :=  if|discn  >  0,2-  fdisc^.oj 


depth  of  hull  axis 
width  of  waterplane 


Q:= 


TA  -  TF 


(*an)2 


disc  :=  ■ 
n  4 


-w 


f  :=  if 
n 


(  2dn  ^ 

dia  >  0 , - ,  1 

n  dia 

V  n  7 


(above,  "I"  is  arbitrary) 
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p  :=  1025  g  :=  9  .81 


Are<x)  :=  linterp(xh,  area,x)  (i  .e . ,  trapezoidal  integration) 


Vol  := 


xhN 


Area(x)  dx 


xho 


MomV  := 


■xhw 


xArea(x)  dx 


xho 


Vol  =  1.2184  x  10 
Zinit :=  P'g  Vol 
h'timt :=  p  g-MomV 
MomV 


LCB  := 


Vol 


=  1.2252  x  108 
M^t  =  1  8985  x  108 


LCB  =  1.5496 


Initial  Estimates: 


Ywp(x)  :=  iintetp(xh , ywp , x)  (implies  trapezoidal  integration) 


These  are 
relative  to 
midships. 


/■xhv. 


Awp  := 


Ywp(x)  dx 


'xho 

TPI  :=  Awp  1025  9.81 
rxhN 


xYwp(x)  dx 


LCF  := 


xho 


Awp 


Awp  =  1483.6074 


TPI  =  1.4918  x  10 


7 


TPI  in  Newtons  per  metre 


lcf  =  6.6629  relative  to  midships. 


IL  :=  (x  -  LCF)  Ywp(x)  dx  II  =  2.0173  x  10 
Jxho 
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MCT  in  Newton-m  per  degree 


MCTd  := 


IL'P 

180 


MCTd  =  3.5403  x  10 


8 


MCTm  := 


iL'P'g 


xhN"xhO 


MCTm=  1.3541  x  108 


MCT  in  Newton-m  per  metre 


Input  Table  for  normal  (down)  force  and  pitching  moment  about  midships 

(since  that  is  where  PC967  calculates  it). 


0 

0 

7.841  -1 05 

1 

-1 .371 4  1 07 

Z  :=  zttig  Z  =  7.8410  x  105 

M  :=  ztrij  M  = -1.3714  x  107 


Mlcf  :=  M  +  Z  LCF 


MLCF  =  -8  4897  x  10 


6 


Estimated  Sinkage  and  Trim: 


wrt  LCF: 

Az:=^ 

TPI 

Az  =  0.05256 

metres 

(since  MCT  is 
based  on  that) 

mLCF 

A9  :=  - 

MCTd 

A9  =  -0.02398 

degrees 

(Trim,  FP  to  AP  in  m) 


Mlcf 

- =  -0.06270  metres,  positive  by  stern 

MCTm 


Iterated  Estimates... 

Definitions:  zT  :=  z^,  +  z 


M7  :=  Mjjjjt  +  M 


ZT  =  1.2330  x  10® 

Mj  =  1.7614  x  108  Relative  to  midships 


TAn(dz,de)  :=  TA  +  dz  +  (LCF  -  xhj 


tan 


dz  is  taken  at  LCF;  also  assume 
small  angle,  etc.  ... 


df^Xjdz.dO)  :=  TAn(dz,de)  -  (X  -  xh^J  |^Q  +  tan^dO  -  y 


discn(x,DA,dz,d9)  :=  - (^(x.dz.d©)^ 
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ywpn(x,DA,dz,d©)  :=  if(discn(x,DA,dz,de)  >  O^-^disc^X.DA.dz.de)  ,o) 


fn(x, DA, dz,d0)  :=  if 


/ 


DA  >  0 ,  - 


2-dn(x,dz,de) 


gn(X,DA,dz,d9)  := 


■, 


DA 


DA 


-,1 


7i  -  I  acos  (fn(x,DA,dz,de))  - 


ywpn(x,DA,dz,de) 


DA 


fn(x,DA,  dz.de) 


arean(x,DA,dz,d@)  :=  if  discn(x,DA,dz,de)  >0,1  g^X.DA.dz.de)  n 

V.  4  J 


2  > 


Dia(X)  :=  linterp(xh,dia,X) 
Solve  Loop:  Given 


and 


&z  :=  Az 


66  :=  A6 


ZT  =  p  g 


My  =  p-g- 


xho 

»xhN 

xho 


arean(x,Dia(x)  ,6z,66)  dx 


x  arean(x,Dia(x)  ,8z,69)  dx 


:=  Find(6z,66) 

Which  gives  the  iterated  estimates . 

Estimated  Sinkage:  Az  =  0.07015 

Estimated  Trim:  ao  =  0.03916 


(may  be  toggled  off) 


metres  (at  LCF) 

degrees 


(These  values  are  for  default  tolerances  less  one  OOM; 
results  using  the  default  tolerances  are  essentially  the 
same  for  these  test  data  .) 

Resulting  draft  at  AP  and  FP: 


AP 

FP 


TAn(Az,A@)  =  10.04089 


m 


TAn(Az,A@)  -  L-|q  +  tan^A6  —  jj  =  8.10851 


m 
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Nomenclature 


Symbols 

a 

Awp 

B{x),B 

c 

CB 

Cf 

cp 

Cp 

Cw 

cz 

CO 

Fh 

Fr 

9 

h 

II 

kw 

K,  M,  N 
L 
r 
s 

S(x),S 

T 

U 

x,  y,  z 
X ,  Y,  Z 
5* 

V 

P 

$ 


radius  of  a  sphere 
waterplane  area 
beam,  maximum  beam 
bottom  clearance 
nondimensional  LCB 
nondimensional  LCF 
pressure  coefficient 
prismatic  coefficient 
waterplane  coefficient 

nondimensional  downforce  based  on  n a2  or  tt r2 
Vermeer’s  sinkage  coefficient,  equation  (13) 
Vermeer’s  trim  coefficient,  equation  (14) 

Froude  number  based  on  depth 

Froude  number  based  on  hull  length 

gravitational  acceleration 

depth  or  height  above  bottom 

waterplane  moment  of  inertia 

waterplane  longitudinal  radius  of  gyration 

roll,  pitch,  and  yaw  moment 

hull  length 

maximum  hull  radius 

squat 

immersed  area,  maximum  area 
draft 

forward  speed 

vehicle  body  axes 

axial,  lateral,  and  normal  force 

boundary  layer  displacement  thickness 

immersed  volume 

density 

roll,  pitch,  and  yaw  angles 
flow  potential 
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Axes  and  Loads 


*  1  ■  ’ 

2,W,Z 


Acronyms,  Abbreviations,  and  Subscripts 

AP  Aft  Perpendicular 

BOW  (subscript)  forward,  generally  FP 

CF  Canadian  Forces 

DRDC  Defence  R&D  Canada 

FP  Forward  Perpendicular 

MID  (subscript)  midships 

LCF  longitudinal  center  of  floatation 

LCB  longitudinal  center  of  buoyancy 
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